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Abstract 

Wc rewrite the Martin-Siggia-Rosc (MSR) formalism for the statis- 
tical dynamics of classical fields in a covariant second order form ap- 
propriate for the statistical dynamics of rclativistic field theory. This 
second order formalism is related to a rotation of Schwinger's closed 
time path (CTP) formalism for quantum dynamics, with the main dif- 
ference being that certain vertices are absent in the classical theory. 
These vertices are higher order in an h expansion. The structure of 
the second order formulation of the Schwingcr Dyson (S-D) equations 
is identical to that of the rotated CTP formalism apart from initial 
conditions on the Green's functions and the absence of these vertices. 
We then discuss self-consistent truncation schemes based on keeping 
certain graphs in the two-particle irreducible effective action made up 
of bare vertices and exact Green's functions. 
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I. INTRODUCTION 

Recently there has been a lot of interest in connecting hot relativistic 
quantum field theory which is important for cosmology and heavy ion physics 
with classical field theory. Although there have been several insightful stud- 
ies in perturbation theory O, a systematic analysis has not been given for 
the Schwinger Dyson equations. One exception is the work of Wetterich 
||, but that work only concerns itself with equal time correlation functions. 
Since formalisms already exist for obtaining Schwinger Dyson (S-D) equa- 
tions for both the quantum [§] and classical problems Q, we thought it 
timely to review the MSR formalism, and cast it in a form so that it could 
be directly compared with the now more familiar CTP formalism Qj. We 
show here that the second order formulation of the S-D Equations in the 
MSR formalism is identical to that of the quantum theory (CTP formalism) 
apart from initial conditions on the Green's functions and the absence of 
those vertices which are higher order in h. 



II. REVIEW OF FIRST ORDER IN TIME MSR FORMALISM 

In the paper of Martin-Siggia-Rose (MSR) Q, an operator formalism 
was developed which allows one to find the generating functional for both 
the correlation and response functions for first order classical field equations 
of the type: 

x{r,t) =A[x(r,t)] , (1) 

where A[x(r, t)] is a local polynomial in the classical field x(r, t). In the work 
of MSR, A[x(r, t)] could contain dissipative terms as well as prescribed noise 
terms. The formalism presented in MSR is first order in time derivatives 
and not apparently covariant. What we will show is that there is a covariant 
subset of the MSR equations in terms of which all the MSR Green's functions 
can be recovered. This subset can be derived from a second order Lagrangian 
formulation which can be related to the h — ► limit of the CTP formalism 
of Schwinger and Keldysh. 

For the statistical classical field evolutions of N interacting classical fields 
(p a , a = 1,2 ... N then x(r, t) is the 2N component field consisting of 4> a and 
the canonical momentum ir a = <j) a . 




If for example we restrict ourselves to cubic interactions, then the vector A 
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is of the form 

Ai = Ci(r,t) + dijXj(r,t) + -eij k Xj(r,t)x k (r,t) , (3) 
where i = 1 . . . 2N. In the MSR formalism one introduces the operator 

x(r, t) 



5x(r, t) 

such that the commutation rule 

[x(r,t),x(r',t)] = S(r - r') , (4) 
is true. Defining an operator Hamiltonian H 

H(t) = J dr' Xi(r' ^A^r' ,t) , (5) 

the equations of motion can be written in the compact form 

x(r,t) = [x(r,t),H(t)} . (6) 
For Eq. (|j) to be true at all times one needs that x satisfies 

dx(r, t) 



dt 



[x(r,t),H(t)] , (7) 



for consistency. Therefore x(r,t) is a functional of x(r, 0) and x(r, 0). The 
formal solution to these equations is given by (in what follows we suppress 
the spatial coordinate r) 

x(t) =U- 1 (t,0)x(0)U(t,0) 

x (t) = f/~ 1 (t,0)x(0)C/(t,0) , (8) 

where 

U~ 1 (t,t ) = Texp[- f H(t')dt']; U(t, to) = T* exp[ ( H(t')dt'\ . 

J to J to 

The meaning of the expectation value (x(t)x(t')) is as follows. Given an 
initial probability function P[x(0)] then 

(x(t)x(t'))= / dx(0) x[t,x(0)] x[t',x(0),x(fi)] P[x(0)] . (9) 
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Wherever x(0) appears, it is replaced by 



i(0) ^ -sjsy ■ 

and it acts on everything to the right of itself. This definition of the ex- 
tended averaging procedure has three important properties. (1) The average 
of a product of x's agrees with the conventional definition. (2)The time 
dependence of (x(t)x(t')) is consistent with the above equations of motion. 
(3) The expectation value of a product of x and x which has an x to the left 
vanishes if P[x(0)] goes to zero fast enough at large \x\. 

The last property is crucial for the tridiagonal form of the Green's func- 
tions and follows from the fact that 

poo J 

/ — (x n P(x)) dx = if lim x n P[x] = . 
J-oc dx \x\-*oa 

Thus in particular in the absence of external sources (x(t')x(t)) = (x(t')x(t')) 
0. The meaning of the hatted operators is understood in terms of the re- 
sponse of the system to an external source. If one changes the Hamiltonian 
by 

SH = J dr'x(r't)Sf(r't) , (10) 
one then finds || that the response of any observable A is given by: 
6{A %sf=o = (T (A(t)x(t>))) 



Sf(t 
Sf(t')6f(t") 



-\ Sf=0 =(T{A(t)x(t')x(t")) , (11) 



where T corresponds to the usual time ordered product operation. Since 
x(t) and x(t') always commute, it is clear that the generating functional for 
both the fluctuation and response functions is given by 

W[f], fj] = In Z[rj, fj] = ln(T(S[rj, f}])) 

S[rj, fj] = exp / dt dr [x(r,t)rj(r,t) + x(r,t)f)(r,t)] . (12) 



In particular the one particle functions are 

8W 5W 

\ v =f,=o = (x(t)) , j-r-r \ v =fj=0 = (»(<)) = . (13) 



5 V (t)>"-"-" — ' 6f)(ty 
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The two particle Green's function matrix is tridiagonal 
G 2 (t) = 

where 



U(t) D R (t) 
D R (-t) 



(14) 



and 



G XX (12) = U(12) = (x(l)x(2)) - (x(l))(x(2)) 



G xx (12) = D R (12) = (T(x(l)x(2))} = @(h - i 2 )(x(l)x(2)} . 



One next doubles the space as in the CTP formalism by introducing the 
field $(r,t) 



x(r, t) 
x(r, t) 



(15) 



Then the commutators of $ satisfy in the( 2x2) space 

[9 a {r,t),9p(r f ,t)] =i{a 2 ) ap 8{r-r') . 
In this larger space the equations of motion are 

$ = [<f>,H] , 

where the "non-Hermitian" operator H has the form 

H = J 7a*a(l)<*l + \llb J d\ $ a (l)*6(l) + ^llbc f *1 $ a $ b ®c , (18) 
and the equations of motion in the presence of an external source J is 



(16) 



(17) 



-ia 2 <t> a = Ja{l) + 7a + lib M 1 ) + ifllbc^c , 



where 



J(r,t) 



rj(r, t) 
fj{r,t) 



(19) 



(20) 



The generating functional W[r], fj] of Eq. |l2| is now written in compact form 
by writing 



5 = exp 



J d r d t H r,t)J(r,t) 



(21) 
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The expectation value of the equation of motion in the presence of an ex- 
ternal source leads to the equation 

- «7 2 Gi(l) = J(l) + 7 2 Gi(l) + ^7 3 [G 2 (H) + Gi(l)Gi(l)] , 

5J(l)5J(2)---5J(m) ' 1 ' 

All the higher n-point functions are obtained by functional differentiation. 
The complete S-D equations are derived in the original MSR paper Q, and a 
path integral representation was derived by Jouvet and Phythian ||. These 
coupled Green's functions are in a 2N x 2N matrix space and obey first order 
in time equations rather than in the usual N xN space obeying second order 
equations. To see how to obtain the second order formalism we will simplify 
the discussion by considering a cubic plus quartic anharmonic oscillator. 



III. SECOND ORDER MSR FORMALISM 

For the anharmonic oscillator with equation of motion 

x + m 2 x + gx 2 + Ax 3 = , 
the above discussion leads to 

H = xp — p(m 2 x + gx 2 + Ax 3 ) . (23) 
This yields the first order equations 

[x,H]=p, ^ = \p,H] = -(m 2 x + gx 2 + Xx 3 ) , 

[x,H] = p(m 2 + 2gx + 3\x 2 ) , ^ = fe, H] = -x . (24) 

dt 

The first order equations for the Green's functions are obtained by tak- 
ing the expectation value of these equations with sources added and then 
using functional differentiation. The expectation value is over an initial dis- 
tribution function /9[xo,Po]- ^ srm pl e example is is to choose our ensemble 
of initial conditions at t = from a thermal distribution of the free Hamil- 
tonian, i.e. 

11 f 

p = Ne~ pH ° , H = -p 2 + -m 2 x 2 , / pdxdp = 1 . (25) 
z Z J 

One can show that out of the twelve nonzero Green's Functions in the first 
order formalism only two are independent. These can be chosen as G xx , and 
Gxp. 



dx 
~dt 

dx 

~dt 



The harmonic oscillator values of these Green's functions are easily de- 
termined from the operator solutions 

x(t) = xq cos mt + — sin mt , pit) = po cos mt sin mt , (26) 

m m 

and the initial density matrix. For the thermal initial conditions described 
above we obtain 

(x(t)x(t')) = cos m(t - t') , (x(t)p(t')) = — sin m(t - t') . (27) 
pm z m 

The response function is 

(T(x(t)p(t'))) = D° ret {t - t') = 9(t - t')— smm(t - t') . (28) 
In general, in the absence of sources 

Gppit -t')=0, Gp X (t - t') = G xp {t' -t) = D° adv (t - if) . 
The independent second order equations (adding sources) are 

d 2 

[-p + m 2 ]x + gx 2 + Xx 3 = j x , 
d 2 

^dt 2 + + ^ gX ^ + 3- ^ x2 ^ = I? ' 

Here j5(0) is treated as an operator -j- when one averages over the initial 
probability function in phase space. 

These equations are derivable from the Lagrangian 



Lmsr = -x I P [-772 + m2 \ x + x ["772 + m2 \ P I + 9% 2 P + ^P ~ Jv x ~ 3xP 



1 ( * , d2 21 ,d 2 
2{ P [ d^ + m] X + X [ dt 2 

(30) 

The vertices of the classical theory are 



"fpxx — ^fxpx — Ixxp — 2(/ , 'Jpxxx — "fxpxx — Ixxpx — ^fxxxp — 6A . 

A formal path integral formalism can be generated for the generating 
functional 

Z[j fj ,j x ] = (j dxdpe-f L ^ dt ) , (31) 
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where the expectation value is over the initial density matrix. 

We want to compare Lmsr with the Lagrangian for the CTP formalism. 
The result of the CTP formalism M is that the action is the difference of 
two terms, one for each branch of the closed time path contour. Explicitly 
for the anharmonic oscillator 



L C tp = ^x + [^ +m 2 ]x + - x_{^ +m 2 ]x. 

+±g[x 3 + - xl] + 1a[4 - xi] - j+x+ + j-x- . (32) 
Introducing the change of variables 

x+ = x + — , x_ = x- — , (33) 
we obtain in this new basis 



r^ 2 ) 1 / if- rfl 

~h — = 2v dt?~^ m xJrX ^^2 +m ^ p)+gxp + Xxp 
gh 2 „ 3 Xh 2 a3 



where 



+ — X P ~ h x ~ 3xP , (34) 



3p — 3+ 3— i 3x — 



2 

This formal manipulation can be justified by first obtaining the SD equations 
in the (+, — ) basis and making a rotation by 7r/4 as discussed in Ref. [1]. 
We notice that there are five extra vertices in the quantum case 

_ gh 2 _ 3A^ 

Ippp — ^ ' Ixppp Tpxpp Tppxp Tpppx ^ ' 

not present in the MSR Lagrangian. Note that h arises from the fact that 
the classical p is the derivative operator -d/dp which for quantum mechan- 
ics becomes jd/dp. It may be noted that apart from these extra vertices in 
the quantum case, h dependence also enters in the initial conditions on the 
Green's functions. This dependence is explicit in the commutator contribu- 
tions and also occers in the h dependence on the initial distribution which is 
now a Bose-Einstein distribution rather than a Maxwell-Boltzmann one. As 
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an illustration consider the quantum version of Eq. (26) for the Harmonic 
oscillator with an initial Bose distribution. Then we obtain 



x Q (t) = x9 cos(mi) + ^- sin(mt). (35) 

m 



From this we find 



7) 

(x(t)x(t')) = (xn) cos(mt) cos(mt') + (— sin(mi) sin(mt') 



where 



and 



-i^sm[m(t-t% (36) 



(4) = ^ + ^bMnPm)-l}-\ (37) 



(p 2 o)=m 2 {x 2 ). 

So we see that h enters in many ways in the initial conditions for the 
Green's functions, and also in the structure of the S-D equations where cer- 
tain vertices are proportional to h?. The way one derives the S-D equations 
from the action is identical for both classical and quantum mechanics. Thus 
we obtain the same structure classically, but there are fewer vertex functions. 

We now derive the S-D equations for a generic cubic self-interacting field 
theory whether classical or quantum. For N interacting scalar fields, we 
introduce the column vector <3? Q composed of = 0, and ^ = 7Tj where 
i = 1,2, ■■■N. We also need the metric g a p which is just a^g as far as 
connecting the ( I >1 and <I> 2 sectors. Then we can write generically for cubic 
interactions: 

C = ^ a (D^) a/3 ^ + i 7 ^ p * a $/V - J a ^ , (38) 

where 

D Q X ap( x ~y)= 9af3p + m 2 ]5(x - y) . 
The generating functional is formally 

Z = (T(exp[J a $ a ])) = e w[J] . (39) 
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Defining the "classical" field 3> a and the connected 2 point function G via 
one has that in the presence of sources 

D„ ^'(12) + -7ofr[ gJ( Y + *"(1)G*'(12) + ^(1)^(12)] 

= «l-2). (40) 

To obtain the S-D equation for the inverse two point function we first use 
the connection between the connected 3-point function and the 1-PI vertex 
function |j| obtained by first Legendre transforming from J to $ and using 
the chain rule 



^d) y d2G(12r ^' (41) 

(r (* a (i)^(2)* p (3))) c = ^M 

= y d4 d5 d6G Qa /(14)G^/(25)G pp /(3,6)r a W(456) . (42) 
where the 1-PI three point function is defined by 

r «* (123) = -^kr ■ (43) 

The usual S-D is obtained by multiplying Eq. ^ on the right by G" 1 to 
obtain 



G-^l, 2) = D-^l, 2) + E a/3 (1, 2) + 7t ^ p ^(l)5(l - 2) , (44) 
and 

Sa/s(l,2) = \j d3d4 lapa GP x (13)G™(U)T(342) Xl/ p . (45) 



We also have 



1^,(1, 2, 3) = laap 5{l - 3)5(1 - 2) + . (46) 
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However S is a proper self energy and can be related Q to the effective 
action l^fG] for the 2-PI irreducible graphs of the theory via 

so that it is just a function of the bare vertices and the full Green's functions. 
Using 

ss{12 ^ - fdjd.5"HT'^ . («) 



(5$ p (3) 7 <5G ctA (5,6) <5$ p (3) ' 
we obtain that 

1^(1,2,3) =7^(1-2)5(1-3) 

-J did5 d6 dl r apa (U5)G pv (46)G aX (57)H(67; 23)^ A; ^ • (49) 
The scattering Kernel H is defined by 

5£ a/3 (12) rr/ 1f)K «x 

^(5,6) = ^ 1256 ^ • 

Self consistent approximations M, || are determined by keeping a certain 
class of graphs in T2 [G] , the sum of all 2-PI graphs made from bare vertices 
and full propagators. For cubic interactions, the Bare Vertex Approximation 
(BVA) is obtained by keeping the graph 

rfW27 4jfc G^(12)G^ m (12))G fc "(12)7,„ m , 

which then leads to the self energy being the one loop diagram, and the 
scattering Kernel being single particle exchange. By using the variational 
definitions of £ and H one is guaranteed an internally consistent approxi- 
mation. As an example of the difference between the quantum and classical 
S-D equations, let us restrict ourselves to a cubic anharmonic oscillator and 
make the bare vertex approximation. The S-D equation for the correlation 
function becomes 

G xx {12) = G° XX (12) + j d3d4Gi° i ) (13)7^G^(34)G fcm (34)7 Zmn G ra (42) , 

(50) 
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which can be identified with Kraichnan's Direct Interaction Approximation 
||. Expanding we can write this in terms of classical and quantum contri- 
butions. We have in symbolic form 

G xx = G° xx + G xx T, c ^Gp X + G x pEftG xx + G x p[Ep + T,p]Gp X . (51) 

The classical contributions to the self-energy in the BVA are 

£g(12) = \^ XX G 2 XX (12) , £g(12) = 7^^(12)^(12) , 
5^(12) =71^.(12)^(12), (52) 
and the quantum contribution is 

S?(12) = ^7^7ppp[G^(12) + G| x .(12)] , (53) 

which has formal order of fo 2 since ^ppp is of order H 2 . 

We have recently shown that the classical BVA gives excellent agreement 
with Monte Carlo simulations in 1 + 1 dimensional cj) A field theory flOR . 

IV. CONCLUSIONS 

We have shown how the classical S-D equations can be derived from the 
quantum S-D equations by comparing the CTP formalism with the MSR 
formalism. Our hope is that this result will allow researchers to make clearer 
the connection between quantum field thory and classical field theory at high 
temperatures. We are in the process of comparing the BVA for the quantum 
and classical cases at high and low temperature to understand for what range 
of temperatures the classical approximation is valid. 
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